Topological order in a correlated Chern insulator 



Joseph Maciejko^ and Andreas Riiegg^ 
^Princeton Center for Theoretical Science, Princeton University, Princeton, New Jersey 08544 j USA 
^Department of Physics, University of California, Berkeley, California 94720, USA 

(Dated: May 16, 2013) 

We study the effect of electron-electron interactions in a spinful Chern insulator. For weak on-site 
repulsive interactions at half-filling, the system is a weakly correlated Chern insulator adiabatically 
connected to the noninteracting ground state, while in the limit of infinitely strong repulsion the 
system is described by an effective spin model recently predicted to exhibit a chiral spin liquid ground 
state. In the regime of large but finite repulsion, we find an exotic gapped phase with characteristics 
partaking of both the noninteracting Chern insulator and the chiral spin liquid. This phase has an 
integer quantized Hall conductivity 2e^/h and quasiparticles with electric charges that are integer 
multiples of the electron charge e, but the ground state on the torus is four- fold degenerate and 
quasiparticles have fractional statistics. 
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Over the past few years, the study of topological in- 
sulators [1, 2] has become one of the most active areas 
of condensed matter research. The prototypical topolog- 
ical insulator is the Chern insulator (CI), first theoreti- 
cally proposed by Haldane [3] , which exhibits the integer 
quantum Hall effect (QHE) in the absence of any exter- 
nal magnetic fields. Many concrete proposals to realize 
this new state of matter have followed Haldane's original 
idea, culminating in the recent theoretical prediction [4] 
and experimental discovery [5] of the CI in thin films 
of Cr-doped (Bi,Sb)2Te3. Although electron-electron in- 
teractions are needed to establish ferromagnetic order in 
these films and break time-reversal symmetry as required 
for the existence of a QHE, this effect is described by 
mean-field theory and the resulting CI state is weakly 
correlated. A burning question in the field of topological 
insulators is whether strong electron correlation effects 
can give rise to qualitative changes in their properties [6]. 

Similar to the spin-polarized v = 1 integer QH state, 
the simplest model of CI such as Haldane's original model 
describes noninteracting spinless electrons and has a Hall 
conductivity Gxy = Ce^/h where the Chern number 
C = ±1. The simplest type of interaction one can add is 
a nearest-neighbor density-density interaction. Many re- 
cent studies have focused on partially filled bands, which 
correspond to a metal in the noninteracting limit and 
lead at certain fillings to a fractional CI [7-11] for strong 
enough interactions and sufficiently flat bands. In this 
paper we consider a completely filled valence band (half- 
filling), corresponding to a CI in the noninteracting limit. 
Analytical and numerical studies [12-14] of the spinless 
interacting problem show that the CI persists up to a crit- 
ical interaction strength beyond which the ground state 
develops charge density wave order. Once the electron 
spin is considered, two possibilities arise. The first is to 
combine two copies of the CI with = 1 for spin up 
and C| = — 1 for spin down. The resulting noninter- 
acting model is the time-reversal invariant quantum spin 



Hall (QSH) insulator [15] in the limit of conserved z com- 
ponent of spin [16, 17]. In this case, the simplest type 
of interaction to consider is an on-site Hubbard interac- 
tion, and the resulting interacting problem at half-filling 
has been the focus of intense recent study [6]. The sec- 
ond possibility is to consider two copies of the CI with 
the same Chern number = = ±1 for both spins. 
The resulting noninteracting model at half-filling is a CI 
with total Hall conductivity cFxy = ±2e^//i. In this paper 
we consider the effect of an on-site Hubbard interaction 
U in such a system. This problem, which has received 
considerably less attention than the interacting QSH in- 
sulator, has been studied recently using the slave-rotor 
technique [18, 19] and the mapping to a spin model [20]. 
In both cases a ground state with topological order cor- 
responding to the Kalmeyer-Laughlin chiral spin liquid 
state [21-23] was found in the phase diagram. This re- 
sult corresponds to the ^ oo limit where charge de- 
grees of freedom are frozen and spin-charge separation 
occurs. Here, motivated by recent work on the inter- 
acting QSH insulator [24, 25] we are interested in the 
possibility of new phases occurring at intermediate U. 
Besides the weakly correlated CI at small U and the chi- 
ral spin liquid at large [/, we find within the Z2 slave- spin 
representation [26-28] a new chiral topological phase at 
intermediate U which we dub the CI* phase by anal- 
ogy with the QSH* phase in Ref. 25. The CI* phase is 
characterized by an emergent deconfined Z2 gauge field, 
a four-fold topological ground state degeneracy on the 
torus, and spin-charge separated excitations with frac- 
tional statistics (Fig. 1 and Table I). 

We consider a class of 2D electron lattice models of the 
form 

H = ^^trr'cl^Cr'a^U^\^^nra-^\ , (1) 

rr' OL r \ OL / 

where (cra) is a creation (annihilation) operator for 
an electron of spin a =t,i at site r, — ^ta^ra is 
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FIG. 1. (Color online.) (a) Topological ground state degener- 
acy in the CI*: energy difference between two different global 
Z2 flux configurations on finite tori with Lx L unit cells. The 
dashed curve is a fit corresponding to an exponential depen- 
dence on L. (b) Spin-charge separation in the CI*: a sin- 
gle hole splits into two anyons, a charge — e particle without 
spin (quasiparticle h in Table I) and a neutral particle with 
spin 1/2 (bound state of quasiparticle h and U). Results are 
obtained from self-consistent calculations on the honeycomb 
lattice model with Hubbard U = 26t and t2 — it. 

the number of electrons of spin a on site r, trr' = tr'r 
a spin-independent hopping amplitude and U is the on- 
site Hubbard interaction. We consider the half-filled case 
(J^^ rira) = and choose the hopping amplitude trr' to 
be that of a spinless two-band CI with (say) Chern num- 
ber 1 in the valence band. Although we will present nu- 
merical results for Haldane's honeycomb lattice model [3] 
later on in the paper, we begin with a more general dis- 
cussion of the possible phases and their expected univer- 
sal properties which are independent of the details of the 
lattice model. 

The model (1) can be first investigated in two extreme 
limits. In the limit = 0, the system is a noninter acting 
CI with total Chern number C = 2 and Hall conductivity 
= 2e^ /h. Because this is a gapped state, we expect 
(Jxy to remain quantized for nonzero but small /7, with 
the ground state evolving adiabatically from a noninter- 
acting CI to a weakly correlated CI. In the limit ^ 00, 
one can derive an effective low-energy Hamiltonian in the 
subspace of one electron per site which, owing to the full 
SU(2) spin rotation symmetry of the Hamiltonian (1), is 
an SU{2) invariant 5* = 1/2 spin model. For a particular 
choice of trr' on the square lattice, such a spin model 
was recently studied numerically [20] and argued to have 
a chiral spin liquid ground state. 

In order to study the possible phases of (1) at interme- 
diate values of f/, we make use of the recently introduced 
Z2 slave-spin theory for correlated electron systems [25- 
28] . This theory is based on the simple observation that 
the Hubbard interaction energy in (1) depends only on 
the total occupation of site r modulo 2^ which can be 
represented by an Ising variable r^, viz. = 1 for un- 
occupied or doubly occupied sites and = — 1 for singly 



occupied sites. The electron operator Cra is written as 

a product of a slave-fermion operator fra and an Ising 
slave-spin r^, 

rot J rot r ' V / 

where and obey the (anti) commutation relations of 
the associated Pauli matrices, so that flips the sign of 
since creating/annihilating an electron changes the oc- 
cupation modulo 2. The Hamiltonian (1) is then written 
in terms of the slave-fermions and slave-spins as 

^ = E E trr'T^r-'fljr'a + | E « + 1) ' (3) 

rr' cx r 

While the original Hamiltonian (1) acts in the Hilbert 
space of physical electron states {IV^phys)}, the Hamilto- 
nian (3) acts in an enlarged Hilbert space which contains 
unphysical states {IV^unphys)} such as states with a single 
slave-fermion /ra/ra = 1 but formally zero/double 
occupancy = 1. When evaluating the partition func- 
tion, one should sum only over physical states. This is 
accomplished by introducing a projector P inside the par- 
tition sum Z = Tr (e~^^P) where P is defined such that 
^IV^phys) = IV^phys) and Pl-^unphys) = 0. Although the 
details of the derivation will be given elsewhere [29], Z 
can be written as the partition function of a Z2 gauge 
theory [30] in 3D Euclidean spacetime with bosonic and 
fermionic matter in the fundamental representation, 

Z= ^e-^[/'/-^-], (4) 

{rnw^J} 

with an action S = Sr^ -\- Sf -\- Sb such that 

— ~ tij fia^ij fja^ 

ij a 

i,j=i — T 

where aij = ±1 is a spacetime Z2 gauge field, z, j are sites 
on a 3D spacetime lattice, and tij is proportional to trr' 
on spatial links and equal to —1 on temporal links. Sb 
is a Berry phase term [31] which corresponds to a back- 
ground static Z2 charge on every site [32]. The constant 
hz depends on the Hubbard interaction strength [/, 

K= ilncoth(^^^ , (5) 

where e is a short-time cutoff [33]. 

We now discuss the possible phases of the gauge the- 
ory (4). In the noninteracting limit [7 = 0, Eq. (5) im- 
plies that ^ 00. In this limit, the gauge fields are 
frozen [29, 34] and can be gauged away, e.g., a^j = 1. 



3 



The only configurations wiiicii contribute to the partition 
function in the K: ^ oo limit are those with the slave- 
spins ferromagnetically ordered, e.g., rf = 1 on every 
site (or gauge equivalent configurations). Equation (2) 
then implies that the electron and slave-fermion opera- 
tors are proportional, and we recover the noninteracting 
CI. In the U ^ oo limit, Eq. (5) implies that ^ 0. In 
this limit the slave-spins can be integrated out. Tracing 
over (Jij on spatial links generates interactions between 
slave- fermions. 



SSf 



rr' r a(3 



frarfr'arfrldTfi 



r' /3t 



where . . . involves higher multiples of four slave-fermion 
operators, while tracing over a^j on temporal links im- 
poses the local constraint frafra = 1 on every site 
r. Therefore the U ^ oo limit corresponds as expected 
to a 5 = 1/2 spin model, which can exhibit a variety 
of ground states depending on the form of trr'- Kef. 20 
argues that a special choice of trr' on the square lattice 
can give rise to a chiral spin liquid. 

We now discuss the regime of large but finite U. In 
general, using the Poisson summation formula we can 
rewrite a gauge theory (here p = 2) as a compact 
U{1) gauge theory coupled to a charge-p integer- valued 
link variable Uij with no dynamics [35], 



-I 



DfiaDfiaDai 



{^i} {^ij} 



where S = Sr- ^ Sf ^ Sn^ Sb with 

ij a 

Sfi — ip ^ ^ TlijCLij , 



where Ciij IS a compact U{1) gauge field and we have 
introduced the external U{1) electromagnetic gauge po- 
tential Aij which couples to the electrons with charge 
e. The summation over riij has the effect of discretiz- 
ing aij in integer multiples of 27r/p. For large but finite 
[/, is small. The slave-spins are gapped and can be 
integrated out perturbatively [34], resulting in a lattice 
version of the Maxwell term for a^j. Because the band 
structure of the slave- fermions is gapped, we expect on 
general grounds [30] that the resulting gauge theory ad- 
mits a gapped deconfined phase for large but finite values 
of U. In the rest of the paper we focus on the physics of 
this phase which we denote the CI* phase. 
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TABLE I. Electromagnetic charge Q, spin Sz, self-statistics 
Oil and nontrivial mutual statistics Ou^ of quasiparticles in 
the CI* phase, where h is the gauge charge vector [36] of 
quasiparticle i. 



In the CI* phase, the gauge field a^j is essentially free 
and we can take the continuum limit a^j a^. Following 
Senthil and Fisher [31], we expect that the Berry phase 
term Sb can affect the position of the phase boundary 
for the deconfined phase but not the universal proper- 
ties of this phase. On the other hand, Sn is essential 
to remember the Z2 nature of the gauge theory in the 
continuum U{1) description. The link variable Uij 
obeys the constraint 9^n^ = to preserve the gauge in- 
variance of the action under a /7(1) gauge transformation 

djyh\ where 



ax, 



of a^. Solving the constraint by 

is a compact U{1) gauge field and passing to a real 
time description, we obtain 

i.e., the (2+l)D level-p BF term [37, 38]. Finally, we 
consider the slave-fermion action Sf. In a hydrody- 
namic approach [36], we introduce a conserved U{1) cur- 
rent = ^e^^^dj^a^ for slave- fermions of each spin 
<j =t, i [39], where ctj^^ctf^ are compact U{1) gauge fields. 
Considering that slave-fermions of each spin form a CI 
with Chern number = = 1, we obtain the effective 
Lagrangian of the CI* phase. 



-CI* 



E 



2tt 



(6) 



which is the main result of this work. From Eq. (6) we 
can extract all the universal topological properties of the 
CI* phase, such as the i^-matrix [36, 40], 
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(7) 



in the basis (a^, 6^, a]^, a^). The ground state degeneracy 
on the torus is | det J^l = 4, as for phases with Z2 topo- 
logical order [41]. Unlike such phases however, the CI* 
phase is a chiral topological phase with Hall conductivity 
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FIG. 2. (Color online), (a) Charge pumping between two 
hexagons threaded by external fluxes and —cj) and sepa- 
rated by D unit cells, (b) Transferred charge as a function of 
(j) in the CI (dashed curve) and the CI* (blue circles), in units 
of the flux quantum 00 = h/e for D = 8. The shaded area de- 
notes the parameter range where a spontaneously created Z2 
vortex pair (semion/antisemion) partly screens the external 
fluxes. 

where t = (0, 0, 1, 1) is the electromagnetic charge vector. 
The i^-matrix (7) can be made diagonal by a GL(4, Z) 
transformation W [36], 

K W^KW = diag(2, -2, 1, 1), (9) 
t^iy^t = (-2,2,l,l), 

which describes the direct sum of the double semion 
model [42] and a CI with Chern number = = 1. 
The quasiparticle excitations in the CI* phase and their 
quantum numbers can be read off from the i^-matrix (Ta- 
ble I). The quasiparticles li (semion), I2 (antisemion), 
and ^3 (semion-antisemion bound state) have the same 
statistics as the excitations of the double semion model, 
but in contrast to their occurrence in spin models, here 
these excitations are electrically charged. Their charge 
and statistics can be understood as follows. An applied 
TT flux in a = 1 integer QH state traps a fermion mode 
with Q = ±e/2, and the resulting flux-charge composite 
behaves as an anyon with statistical angle f [43]. The 
semions/antisemions in our model likewise arise from in- 
gap modes of the slave-fermions with total Chern num- 
ber 2, which are bound to dynamical Z2 vortices. The 
quasiparticles I4 and ^5 are identified with the original 
electrons, and are topologically equivalent to the vac- 
uum [44]. Although the quasiparticle spin Sz in Table I 
is only meaningful in the presence of U{1) spin rotation 
symmetry, the other properties of the CI* phase are ro- 
bust under perturbations which break this symmetry but 
do not close the bulk gap. 

The existence of protected edge excitations can also be 
inferred from the i^-matrix [36]. K has an imbalance of 
positive and negative eigenvalues [Eq. (9)], which trans- 
lates into an imbalance of right-moving and left-moving 



edge excitations on a manifold with boundary. Such a 
chiral edge cannot be completely gapped out by local 
perturbations [45]. However, because the slave-spins are 
disordered in the CI* phase ((rf ) = 0), the electron and 
slave-fermion operators are not proportional. We there- 
fore expect that the electron spectral function on the edge 
is gapped while correlation functions of on-site electron 
bilinears cl^Cr/s or CraCr/3 are gapless, as in the QSH* 
phase [25]. 

To further study the CI* phase we use the self- 
consistent mean-field framework in the Hamiltonian lan- 
guage [25, 27]. In this scheme, the bond energies of the 
free slave-fermion model provide the exchange couplings 
Jrr' = ^ai^rr'ifrafr'a) + c.c.) of the transversc-field 
Ising model, whereas the correlation functions grr' = 
(r^T^/) of the transverse- field Ising model renormalize 
the hopping parameters of the slave-fermions [46]. Per- 
forming unrestricted mean-field calculations in the CI* 
phase on finite systems, we find a four- fold topological 
ground state degeneracy on the torus [Fig. 1(a)] as well 
as spin-charge separated excitations above the ground 
state [Fig. 1(b)]. For these numerical results, we studied 
the spinful Haldane-Hubbard model on the honeycomb 
lattice with U = 26t and a complex second-neighbor hop- 
ping t2 = it where t is the nearest-neighbor hopping. 

Although the Hall conductivity (8) of the CI* is the 
same as in the weakly interacting CI, interesting differ- 
ences appear beyond the linear response regime (Fig. 2). 
We study the induced charge in a charge-pumping setup 
where an external flux is threaded through one hexagon 
and a flux —0 through an other hexagon separated by 8 
unit cells. While for small fluxes the pumped charge fol- 
lows the expected relation dQ = C0/0o with C = 2, we 
observe a sudden jump of AQ = e for slightly larger 
then 0o/4. This jump is associated with the spontaneous 
creation of a Z2 vortex pair (semion/antisemion) which 
is localized at the fiux-pierced hexagons and partially 
screens the external fluxes. A second jump of AQ = e 
occurs at a flux slightly less than 30o/4. This behavior 
is in sharp contrast to the charge pumping observed in 
the conventional CI (dashed curve) where the screening 
of the fluxes by Z2 vortices cannot occur. 

In summary, we predict that a new chiral topological 
phase of matter, the CI*, can arise when a spinful CI is 
subjected to a strong on-site Hubbard interaction. Some 
properties of the CI* are reminiscent of noninteracting 
topological phases (integer quantized Hall conductivity 
and integer charge of quasiparticles) while others indi- 
cate the presence of topological order (nonzero ground 
state degeneracy on the torus and fractional statistics of 
quasiparticles). To test these predictions, a determina- 
tion of the phase diagram of Hamiltonian (1) via exact 
numerical diagonalization or density-matrix renormaliza- 
tion group studies is desirable. The interesting question 
of whether the chiral spin liquid of Ref. 20 can be reached 
via a continuous transition from the CI* can perhaps be 
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studied analytically [47] using the topological field theory 
(6) as a starting point, and is left for future work. 

After the completion of this work, we learned of a re- 
lated study of the CI* phase by Zhong et al. [48]. 
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